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Abstract 

We introduce a novel method for the renormalization of the Hamil- 
tonian operator in Quantum Field Theory in the spirit of the Wilson 
renormalization group. By a series of unitary transformations that suc- 
cessively decouples the high-frequency degrees of freedom and partially 
diagonalizes the high-energy part, we obtain the effective Hamiltonian for 
the low energy degrees of freedom. We successfully apply this technique 
to compute the 2-loop renormalized Hamiltonian in scalar X 4> 4 theory. 
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1 Introduction 



The subject of renormalization for the Hamiltonian operator has been looked 
into extensively for the past decade [JD. Yet it is hardly being used in any 
practical computation in field theory, since usually all the intermediate steps 
are non-covariant and, generally, rather complicated compared to the standard 
approach. Nevertheless the concept of "integrating out" of the high frequency 
degrees of freedom has much more physical meaning precisely in the Hamiltonian 
framework. One way to present it would be to set up the whole procedure as a 
Born-Oppenheimer approximation used in atomic physics H. Another way to 
look at it would be to introduce a unitary transformation in order to decouple 
"high" and "low" modes and then look at the low-energy part of the spectrum 
H . In this letter we want to introduce a renormalization technique in the spirit 
of the Wilsonian renormalization [|J, i.e., integration of fast degrees of freedom, 
appropriate for the Hamiltonian (Schrodinger) formalism ||. It is similar in 
essence to both || and || , but appears to be very different in practice. Basically 
what we want to see is how the high frequency degrees of freedom modify the 
low energy Hamiltonian operator, or equivalently, what Hamiltonian in term 
of low frequency modes produces the same low energy physics as the original 
Hamiltonian. This technique was successfully applied in a previous paper by 
the authors to the 1-loop renormalization of Abelian and non-Abelian gauges 
theories ||. 

For the sake of explicitness let us consider a system described by a Hamil- 
tonian H(A) defined up to the momentum scale (cut-off) A; that is, H(A) 
produces finite results. Now if fj, is a lower momentum scale, we want to find 
a new Hamiltonian operator H(p) that generates the same results for all the 
physical processes which do not involve momenta greater than \x. The definition 
of H(fi) is simple: is the projection of the original Hamiltonian H(A) onto 

the low frequency subspace, i.e., the high frequency vacuum, 

H(n)=P low H{A)V low . (1) 

The claim is that, at least perturbatively, it is possible to decouple the low fre- 
quency modes from the high frequency modes and thus give sense to the notion 
of "low frequency subspace" . In essence we will show that after a suitable uni- 
tary transformation we can partially diagonalize the Hamiltonian and construct 
the vacuum for high momentum modes. Therefore we can rewrite eq. (Q) as 

H{n) = (0 M gh\U\A,n) H(A) U{A,fx)\0 Mgh ), (2) 

Now let us explain how to construct the unitary operator U(A,/i). As we 
have explained before, the main goal is to identify the ground state of the high 
frequency modes and then project the whole Hamiltonian onto this state. So if 
we can find a unitary transformation that: i) separates the low energy modes 
from the high energy modes and, ii) diagonalizes the high energy part of the 
Hamiltonian (for example in terms of the creation and annihilation operators 
of the high frequencies), we are done: the ground state will be then the state 
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annihilated by all the high frequency annihilation operators. This task is of 
course difficult, but in fact we need less than that. Since we only need to identify 
the vacuum state of the high momenta Hamiltonian and not the whole spectrum, 
only a partial diagonalization is enough. In fact, only those non-diagonal terms 
containing purely creation operators or purely annihilation operators have to 
be removed from the Hamiltonian. Other non-diagonal terms containing both 
creation and annihilation operators do not change the vacuum state of the theory 
(once the aforementioned terms have been removed). 

We now show how this work in more detail. Suppose that after a unitary 
transformation we bring a Hamiltonian to the form Hdiag + V where V contains 
only terms with at least one creation operator and one annihilation operator. 
That is, V can be written as V — a\Mija,j where M is an arbitrary operator 
and i,j are generic indices. Then, standard perturbation theory tell us that the 
correction to an arbitrary state \n) is given by 



Therefore, if \n) = |0), the vacuum state, it is annihilated by V and there is 
no correction to it at any order in perturbation theory. The ground state is 
unaffected by V. 

Now we can proceed to find the unitary transformation of eq. (|^) . First we 
split the original Hamiltonian into four pieces, 



Here H± contains only the modes with momenta less than /i, H 2 is the free 
part for all the modes with momenta greater than \i and Va + Vb contains 
mixing terms and all non-diagonal high-momentum operators; Va contains the 
"pure" terms that have only high frequency creation operators or high frequency 
annihilation operators, but not both, and Vb the "impure" remaining terms 
(we assume here that Va and Vb are normal-ordered with respect to the free 
perturbative vacuum). Then we break the unitary operator U in a product 
series, U — Uo U\ ■ ■ ■ U n ■ ■ ■ and we compute all the terms successively. The 
objective of each individual U n is to partially diagonalize the Hamiltonian, at a 
given order in A (a generic coupling of the theory) and /u/A. Let us show how 
to achieve this. 

We proceed iteratively by first diagonalizing the Hamiltonian at leading order 
in A, up to the desired accuracy in fi/A. Consider first a unitary operator U 
written as e lS1 ° . Now we perform a unitary transformation on equation (Q) , 
expanding in powers of O (in the general case f2 is at least of order A so at a 
given order only a finite number of terms are needed): 



e~ in ° {H x +H 2 +V A + V B )e iQo = H Y + H 2 + V A + V B + i [H u Qq] + 



As we have explained above, we want to eliminate the "pure" mixing terms Va, 
while we generate an expansion in /i/A. Thus, we impose the following condition 




(3) 



H = + H 2 + Va + V B - 



(4) 



+ i [H 2 ,n ] + i [V A , Slo] + i [Vb, O ] • • • (5) 
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on fio, 

i[H 2 ,n } + V A = 0. (6) 

This equation can be solved perturbatively and since commutators with 
H 2 generate time derivatives of the high frequency fields we have the desired 
expansion. Notice that equation (^) defines flo up to terms that commute 
with H 2 . This ambiguity corresponds to the redefinition of the low energy 
Hamiltonian Hi by a unitary transformation. We will assume some kind of a 
"minimal" scheme - namely that fio does not have a part that commutes with 
Hi. 

After Slo is chosen to cancel Va in the effective Hamiltonian, a new mixing 
term of order A appears from i [Hi, f2o]. However this new term is of a higher 
order in fi/A and will be eliminated by the next unitary transformation Ui = 
e tQl . Explicitly, 

e^e-^iHi +H 2 + V A + V B )e in °e mi = Hi+H 2 + V B + 

i [Hi,Qq] + i [Hi, fii] + i [Hi, Hi] + second order terms • • • (7) 

and we now choose f2i so that 

%[Hi,n ]+i[H 2 ,n 1 ] = o. (8) 

Using equations (0) and (||) we obtain: 

e- mi e- m °(Hi +H 2 + Vi 2 )e in °e mi = Hi+H 2 + V B + ^ [V A ,n ] + 

i [v B , n \ + 1 [v B , nj + i [Hi, ni] - i [[Hi,n ] , o ] - 

- [[H u fi ] , Hi] - g [[Hutoi] , Hi] + higher order terms • • • (9) 

Now it is easy to deduce the logic of the procedure: the next step is to 
introduce 2 in order to cancel i [Hi, Hi], which is of higher order in /x/A, and 
continue with the same process until we have attained the desired accuracy. 
Notice that by construction 

[H 2 ,Q n+ i] = -[Hi,Q n ] (10) 

and Hi w A, Hi w /i, then 

^n+l ~ ^n, (11) 

and any new il is smaller than the previous one by a factor of /x/A. 

Of course so far we have only eliminated the high momenta degrees of free- 
dom up to the first order in the coupling constant. Requiring the absence of 
the A 2 - order mixing terms will lead to the introduction of a whole new series 
of unitary transformations, and the same arguments can be applied to them. 

As a final remark, notice that the "partial" diagonalization at first order in 
A defines a correct low-energy effective Hamiltonian (after projecting onto the 
high frequency vacuum) which is valid up to the order A 3 . The reason is that 
since all the relevant terms of order A have been cancelled, the introduction of 
the new series of O's of order A 2 will only modify the Hamiltonian at order A 4 . 
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2 A quantum-mechanical example 

In this section we illustrate the ideas of the previous section by studying in 
some detail a quantum mechanical "toy model" where the main ingredients of 
our technique can be presented in a simpler context. 

Consider the following Hamiltonian of two coupled anharmonic oscillators, 

H = y (P 2 + Q 2 ) + y (P 2 + q 2 )+X(Q + qf (12) 

where lo\ << oj2- Now we want to find the effective Hamiltonian for P and Q 
resulting from "integrating out" the high frequency modes p and q. We choose 
to compute the effective Hamiltonian up to to order A 3 and (wi/u^) 3 . 

According to the discussion of the previous section we divide the Hamiltonian 
in four pieces as in eq. (f|), 

H = Hi+H 2 +V A +V B (13) 
where the various terms have the form: 

Hi = f (P 2 + Q 2 )+A^ + 3Q 2 + Q 4 ), 
H 2 = uJ2 a a, 

V A = A {|« 2 + | at2 + i fl4 + i« t4 + V2Q(3a + 3at+a 3 + a t3) + 
3Q 2 (a 2 + a t 2 )+2lQ 3 (a + at)}, 

V B = A{3a t a + a t a 3 + 6a t2 a 2 + a t3 a + 3\/2Q(a t a 2 + a t2 a) + 

6 Q 2 ata} . (14) 

In eq. ( |l4| ) , a and a) are the annihilation and the creation operators respectively 
of the high frequency states, defined through a = (p + iq)/y/2 and a 1 = (p — 
iq)/V2. 

Now we are ready to perform the first unitary transformation Uq — e lfl ° . As 
explained in Section 1, f2o is chosen in order to satisfy equation (Q). Using the 
expressions given in eq. ( ft4| ) we obtain: 

O = i\—l-(aV-a 2 ) + V2Q(3{a^-a) + -{a^-a 3 )) + 

\ Q 2 (at 2 - a 2 ) + 2§Q 3 (at - a) + ^ - a 4 )| . (15) 

After performing the unitary transformation with the above expression for 
fio, we obtain the leading order approximation of the diagonalization of the 
whole Hamiltonian: we have eliminated the non-diagonal "pure" terms Va at 
the expense of creating new ones in [Hi , CIq] . However those new terms are of 
order of uj\/uj2 as can be seen clearly from eqs. (|l4|) and (|l5|). 
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The following step is then to cancel this newly created non-diagonal terms 
with the next unitary transformation U\ = e iSl1 where fix satisfies equation (^). 
A straightforward calculation gives 

fli = i\^(^^ 2 +a 2 )+tV2p(^{^ + a) + ^ 3 + a 3 ^j + 

i\QP{a T2 + a 2 ) -i6V2Q 2 P(a^ +a)J. (16) 

Finally, to reach the desired accuracy in — we need also n 2 , defined through 
the relation [HI, Q a ] + [H 2 , 2 ] = 0. We find 

n 2 = -A^ jl2(a t -a)-i\/2Q^3(a t -a)-^(a t3 -a 3 )^- 
i^j (P 2 + Q 2 ) (a) 2 -a 2 )+i 12%/2 QP 2 (a t - a) i &V2 Q 3 (a* - a) j . (17) 

Now we are ready to compute the "low energy" effective Hamiltonian. Re- 
calling the discussion of the previous section, we have the following expression 
for the effective Hamiltonian, up to the order A 3 and l/wf : 

H ef f = (o| jifi + h 2 + v b + l -[v A , n ] + i[v B ,n ]+ i\y B ,sii]+i\y B , n 2 ]- 

i[[Hi,n ],no] - ^[[ifi^o],^] - i[[T^,Oo],f2o] - ~[[Vk,fio],«i] - 
i [[Vb , no] , n ] -Wb, n ] , ni] - ^ [[[ifj , n ] , n ] , n ] | |o). (is) 

Here |0) is the vacuum associated to the operators a and leaving the P and 
Q operators untouched. For the sake of completeness we have written the whole 
expression derived from the rules described in Section 1. However, some of the 
terms vanish upon projection. 

Finally after a long but straightforward computation we have the desired 
low frequency Hamiltonian : 

u\ , o2 , rV\ , x {3 2 4 \ 2 ( 211 , 29wi 



153 lu 2 ( 1 ui ^.wfX^o 166 w? o 

1 + -31— + 45— + 144^ Q 2 - — 1 P 2 + 

4 \ w 2 w2 2 / 9 

162 i ^ Q P - 8l4 Q 2 P 2 + 288 z ^ Q 3 ^ + f 36^i - — 

0^2 1^2 W 2 V w 2 w 2 

72 4 <^ _ it cA + * (*»_^ _ m ^ + (™± 

2 C^/ V W 2 W 2/ V W 2 W 2/ 

^| g 8 ) + 0(a 4 ,i/o4). (19) 
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The skeptical reader can verify, using standard Rayleigh-Schrodinger pertur- 
bation theory, that the spectrum of the effective Hamiltonian ( [l9| ) is the same 
as the low energy spectrum (i.e., the part of the spectrum that remains finite 
if u>2 — > oo) of the original Hamiltonian (O). 



3 Scalar Field Theory 

In this section we apply the same formalism to the case of a scalar field with 
the quartic self-interaction, A </> 4 . We will determine the effective Hamiltonian 
H up to the two-loop order. 

Before proceeding to the loop calculation we have to explain how renormal- 
ization is performed in our formalism. As was shown in Section 1, by projecting 
to the vacuum state for the "high" momentum modes we obtain the effective 
Hamiltonian H e ff. Since our system is supposed to have had an ultraviolet 
cut-off A from the very beginning, H e jf will explicitly depend on this UV scale. 
The renormalization procedure consists of modifying the original Hamiltonian 
H by introducing renormalization Z-factors that depend on the UV cut-off A 
and some arbitrary "renormalization scale" M. Each of the Z's depends on A 
and M in such a way that the effective Hamiltonian obtained from it does not 
depend on A; in fact it has to look exactly like the original one except that all 
A's should be replaced by \x - the scale down to which we are integrating out. 

Let us start from the "bare" Hamiltonian 

H = Jd 3 x Qn|(z) + i$(a;) (-v 2 + m 2 ) $(*) + A$ 4 (z)) , (20) 
and introduce a Z factor for each composite operator, 

(21) 

Each of the Z-factors has a perturbative expansion in A, where we are assuming 
that A has been defined by some renormalization prescription at the renormal- 
ization scale M (A is the "renormalized" coupling in the language of the standard 
renormalization group). Generically: 

Z = 1 + /x(A)A + / 2 (A)A 2 + / 3 (A)A 3 + • • • . (22) 

The functions f n will be chosen order-by-order from the requirement that after 
integration of the modes from /i to A all the corrections sum up in such a 
way that Z(A) — > Z(fi). When doing the one- loop corrections one can therefore 
assume that all the Z's are initially 1 and choose the corresponding f's from the 
condition that high-cutoff dependence be cancelled after computing the He// 
to one loop. For the second loop we use the one-loop Z-corrected Hamiltonian 
and determine /„ to the next order in A by the same procedure. 

According to our philosophy we have to identify the purely "low" part Hi, 
the free "high" momentum part H2 and the interaction terms Va and Vb ■ First 



x 
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we split the original field into its low and high frequency components: 

ct>{x) = ^$ fc e lte , X {x)= ]T ® kelkX - ( 23 ) 



k<fj, fi<k<A 



By virtue of the equations (p3[), J (f>(x)x{x) = and the original Hamiltonian 
( p0| ) can be rewritten as sum on the four pieces H = Hi + Hi + Va + Vb'- 

Hi = |Q7r2(x) + i0(x)(-v 2 + m 2 )0(x)+A0 4 (^ , (24) 
Hi = J Y^(*) + ^)(-V 2 W)xOr)), (25) 

Va + Vb = Ay"( x 4 (x)+40x 3 + 60V+40 3 x). (26) 

In order to separate Va and Vb, according to the discussion of Section 1, we 
introduce the second-quantized form of the high-momentum modes %, 

X (x)= J2 ~iw( akelkx+a l e ~ lkx ) ( 2? ) 

fj,<k<A * k 



Here oj k = Vk 2 + m 2 and can be taken ~ |fc| for /i >> m. The operators 
a k and a£, satisfy the standard commutation relations afc,at] = Skp- Upon 
substituting this definition in the expressions (Bq) and (26) and normal-ordering 



the result with respect to "high" modes we can identify all the terms that contain 
only creation or only annihilation operators; these terms form Va- The rest of 
the interaction part forms Vb- 

H 2 = ^ u ka\a k , (28) 

/i<fe<A 



= v> e^ k +P+"+ r ^a k a p a q a r | 2 ^x)e< k +P+^ x a k a p a q | 

4^/U} k UlpUl q UJ r yjl LU k UJpLO q 



V B 



3 2 (x) e^ fc + y > x a fe a p 4 3 (a;)e 4fc:l; a fc 3 1 a fc a_ fc 




h. c. , (29) 



y/U>kUJpUJ q LJ r ^/U> k UJpUJ q U) r 

xytP+i-V* a\a p a q \ 



UJkUpUlq 



(30) 



Notice that already at this stage we have some contributions to the effective 
Hamiltonian due to the normal-ordering of the terms of type 4> 2 (x)a k aj 1 ; this 
term explicitly depends on the UV cut-off A and can be included into Hi , 

5Hi=6\ J2 />(*)« 3 A^=^ />(*), (31) 

H<k<A L ° k J J 
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where we have made the standard replacement ^2 k — > J mjw- Equation ( |3l| ) 
is, of course, the standard "tadpole" , one-loop mass renormalization. 

In order to pick up the only other one-loop contribution, the coupling con- 
stant renormalization, we have to follow the general procedure of Section I and 
determine Oq f° r the first unitary transformation. Using equation (M) we deduce 



o = H)a]T 



4 ct) 3 (x)e lkx a k + 3 4> 2 {x)e i{k+p ^ x a k a p + 3 1 a k d- k 

^/20J^ LU k ^/UJ k U) p (LU k +u) p ) Auj p u) 2 k 



2<f>(xy ( - k +P+^a k a p a q + e l( - k +P+^ x a k a p a q a r _ 1 ^ 



y/2 UJ k UJ p UJ q (uJ k + UJ p + UJ q ) 4 ^/LO k UJ p UJ q UJ r (uJ k + L0 p + L0 q + LO r ) 

To find the correction to the Hamiltonian we expand the unitary transfor- 
mations and project onto the high energy vacuum. The contributions to H e ff 
are exactly the same as in equation (|l~8l). 

Now we have to determine the potentially divergent contributions (if A — > oo) 
that will emerge after projection to the vacuum state. In doing so the following 
naive power-counting rule is useful: 

a) Each Q n brings fc( n+1 ) (momentum) to the denominator. 

b) Each contraction contributes with k^ 1 . 

c) For each loop integration we have a contribution of fc 3 . 

Therefore, each term for which the total degree of divergence of a term is 
greater than or equal to zero, has to be taken into account. Of course, there 
may be overlapping divergences in computing the proper expressions, but those 
have to be studied individually, term by term. 

By inspection of equation (|l8|), is easy to see that at this order there are 
no contributions to the terms II 2 and $(— y 2 )$ (i. e., = = 1) and 
the only other one-loop term is the coupling constant renormalization. This 
contribution comes from the commutator [V^4,Oq] when contracting twice the 
high energy fields. The rest of the terms are not important at this stage: the 
commutators [Vg,f2o] an( i [Hi,£Iq] are zero when projected onto the vacuum. 
After projecting onto the high frequency vacuum we have 

dp dq 1 



(33) 



(27T) b LU k Up{UJ k + Up) 



9A 2 , /A 
tt hi — 

2tt 2 VM 



cj> 4 (x) d 3 x 



Expressions (|3lj) and ( |33| ) give the two contributions to H e ff at one-loop. From 
them we can deduce the renormalization factors Z m and Zy. 



(A 2 -M 2 ) 9A 2 / A i 



This finishes the one-loop renormalization procedure. (In refs. [g, [7j there 
is an extra one-loop counterterm of order A due to the definition of the "effec- 
tive" wave functional in the Schrodinger representation). Before going to the 
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next order in coupling, we should notice that there is a two-loop contribution 
coming from the same term [Va, S^o] m (P^|), when contracting the high energy 
fields three times. Its leading divergence is quadratic and it has a logarithmic 
subleading divergent part that gives rise to the two-loop wave-function renor- 
malization, 

- -^\[ d 3 r<t,(r)<t>(-r)\ [ p** , 

V (27Tj y [J J J WpW g UV-p-g(Wp+Wg + i*V-p-g) 

= -^(2 1n2-l)A 2 |i0 2 (,)d 3 a; + ^ln^ j \(y^)f (35) 

Extra contributions to H e f f at two- loops come from the terms 

5H=(0\ f-l[[H 1 ,n ] ) il 1 }+i[V B ,Cl 1 }-^[[H 1 ,ft 1 },n 1 ]\ |0). (36) 

However using the power counting rules described above we can see that the 
only divergent contribution comes from the first term. Following equation (0) 
we determine the operator Oi : 



Oi = {-i)X 



'4 [H^^ix)] e ikx a k 3 [H u <t> 2 {x)] e^ k +P>a k a p 



/2uj k Lu k ^2u>kUJ p (uJk + Up) 



2 



2[H 1 ,^x)}e^ k +P+^a k a p a q | P+«+ r > a k a p a q a r ^ 

yJ2u) k U) p 0J q {u) k +LO p +UJ q ) 2 A^jLO k UJ p UJ q UJ r {uj k + UJ p + UJ q + UJ r ) 2 

and after evaluating a momentum integral similar to the one in eq. 
finally get 

*W$ = 4^(±)l[±4*x\ (38) 




M/ [J 2 

It is important to point out that the corrections to {\/4>) 2 , eq. (|35|), and 7r 2 
eq. ([38|), have equal magnitude and opposite sign as it must be, due to the 
Lorentz covariance of the theory. In fact, since is represented by j%, then 
if (j) 2 gets corrected by Z then n 2 ^ should change by ^, so that the equal time 
commutator is preserved. Now we can define the two-loop wave-function renor- 
malization factors and Z v as 

Z t = 1-f (39) 

Z^ 8 7T 4 \ fl 



Z* = _ = l + -- s ]n(-). (40) 



Finally there is one more two-loop contribution that comes at the order A 3 
and renormalizes the </> 4 (cc) term. It comes from the terms 

SH = (0| (~ [[V A , Oo] , Oo] - i [[Vs, O ] , fio]) |0) (41) 
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After a somewhat tedious computation we get 



SB 



(2-loop) 



27 1 



X 6 i In 



In 



^(x) d 3 x. (42) 



To compute the entire two-loop correction we have to add to this result the 
contribution of the one-loop counterterms ,i.e., the one-loop terms with the 
order A contributions to the Z-factors. Thus, we can finally deduce the value of 
the Z\ factor at two-loops, 



1 



2tt 2 



log(A/M) A 2 



(— 

Wtt 4 



(log(A/M)) 2 



51 

4^4 



log(A/M)) A 3 . (43) 



From equations 
theory, 

/3(A) 



I) and (f43|) we get the correct two-loop /3-function of the 



dX 



A 



dlogM 
2tt 2 4tt 4 



((9a(AZ a )| a ^) _1 A 



dZ x 



logM 



A. A 



and the correct two-loop anomalous dimension 7: 



laiogZ, 



2 dlogM 



16tt 4 



-A 



(44) 



(45) 



In summary, we have described in this letter a novel perturbative technique 
of renormalization in the Hamiltonian (Schrodinger) formalism. We have showed 
that this method successfully gives the two-loop renormalized Hamiltonian for 
a scalar Field Theory with a quartic potential. Moreover, in a previous work || 
we have shown that this technique also gives a consistent one-loop renormaliza- 
tion for the Quantum Electrodynamics and the Yang-Mills Hamiltonians where 
results similar to those from Coulomb gauge covariant calculations were found. 
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